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^ • Abstract 

O^ 

^H '. A new model to describe the dynamics of particles undergoing diffusion in general 

^^ relativity is proposed. The evolution of the particle system is described by a Fokker- 

Planck equation without friction on the tangent bundle of spacetime. It is shown 
\^ . that the energy-momentum tensor for this matter model is not divergence-free, which 

^T) \ makes it inconsistent to couple the Fokker-Planck equation to the Einstein equations. 

^" ■ This problem can be solved by postulating the existence of additional matter fields in 

spacetime or by modifying the Einstein equations. The case of a cosmological scalar 
. field term added to the left hand side of the Einstein equations is studied in some details. 

f— ^ ■ For the simplest cosmological model, namely the flat Robertson- Walker spacetime, it is 

shown that, depending on the initial value of the cosmological scalar field, which can be 
identified with the present observed value of the cosmological constant, either unlimited 
expansion or the formation of a singularity in finite time will occur in the future. Future 
^ ■ collapse into a singularity also takes place for a suitable small but positive present value 

H , of the cosmological constant, in contrast to the standard diffusion-free scenario. 

1 Introduction 

Diffusion is one of the most fundamental macroscopic forces in nature. It is the driv- 
ing mechanism of many dynamical processes in physics, such as heat conduction and 
Brownian motion, and it is often invoked to describe transport phenomena in biology 
or even in social sciences. Given the wide variety of phenomena that display diffusive 
behavior, it is somehow surprising that a consistent theory of diffusion in general rela- 
tivity is at present still missing. Besides the evident theoretical physics motivation for 
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developing such a theory, diffusion may also play an important role in the large scale 
evolution of matter in the universe. 

The relativistic theory of diffusion processes has received considerable attention in 
recent yearqj. Most works so far have been concerned with the justification and the 
analysis of the stochastic differential equations that describe diffusion at the micro- 
scopic level. The theory started with the pioneering work of Dudley [11], who showed 
that Lorentz invariant Markov processes do not exist in Minkowski space, while only 
one such process with continuous paths (Brownian motion) can be defined on the tan- 
gent bundle (the relativistic phase space). Dudley's process was generalized to curved 
manifolds in [17] and studied on some exact manifold solutions of the Einstein equa- 
tions of general relativity in [H [5l [T8l [20l [22] . These works however are only concerned 
with describing the diffusion of particles on a given background spacetime, and no at- 
tempt is made to relate the curvature of spacetime with the dynamics of the matter 
undergoing diffusion. 

A different approach to diffusion is based on the kinetic Fokker-Planck equation, 
which is the partial differential equation satisfied by the (one-particle) distribution 
function /. Although the Fokker-Planck equation is formally equivalent (via Ito's 
formula) to a system of stochastic differential equations, it is often easier to justify the 
evolution of the function / than the dynamics of the microscopic stochastic process 
(for instance by an approximation procedure starting from a master or Boltzmann 
equation |24[ I27j). This point of view was taken in [1], where a relativistic Fokker- 
Planck equation previously derived in [12] by stochastic calculus methods, was re- 
discovered as the most natural second order transport-diffusive equation that is Lorentz 
invariant in the absence of frictioro. Moreover, the associated stochastic differential 
equations, when expressed in terms of the proper time, reduce to Dudley's process |21j . 

When considering diffusion in general relativity, the need to rely on the macroscopic 
continuum description provided by the Fokker-Planck equation becomes transparent. 
In general relativity, in fact, the geometry of spacetime is not given in advance but it 
is determined by the Einstein equations, 

Rt.u-^gt.uR = Tf,„ {8ttG = c=1), (1) 

where the matter fields enter through the energy- momentum tensor Tf^u- In kinetic 
theory the energy-momentum tensor is given by a suitable integral of the particle 
distribution /, which in the presence of diffusion is the solution of the Fokker-Planck 
equation. In its turn, the metric, i.e., the solution of the Einstein equations, appears 
in the Fokker-Planck equation. Thus the consistent description of matter undergoing 
diffusion in general relativity requires to study the coupled system of Einstein and 
Fokker-Planck equations. 

However it is not difficult to see that a matter distribution undergoing diffusion 
cannot appear as the only source in the Einstein equations ([1]). The reason is that under 



^See [T], [3]-[B], [5]-[Il], [II]-[11]- The review [12] contains more references and an historical introduction 
to the relativistic theory of diffusion. 

^ There exist other models in the literatm^e that are named "relativistic Fokker-Planck equations" [^112). 
but none of them is Lorentz invariant in the absence of friction. 



the action of the diffusion forces the kinetic energy of the particles is not preserved, 
and accordingly the energy-momentum tensor T^^'^ does not fulfill the compatibility 
condition V ^T^'^ = required by the Bianchi identities and the Einstein equations. A 
possible solution to this problem is to assume the existence of additional matter fields 
in spacetime which exchange energy with the particles undergoing diffusion. Another 
solution, which is studied in some details in the present paper, is to add a cosmological 
scalar field term in the left hand side of ([T|), leading to the following modification of 
Einstein's equations 

R^lv - -^OfiuR + ^'Qfj.u = T^u- (2) 

The evolution of the scalar field (p is determined by the particle distribution function 
/ through the equation V^cj) = V^T^^, which follows by ([2]) and the Bianchi identity 
V^{Rfjiy — -^QfiuR) = 0. In fact, it will be shown that the energy-momentum tensor 
associated to a solution of the Fokker-Planck equation satisfies 

V^T^"' = 3a r, 

where a > is the diffusion constant and J^ is the current density of the matter. Thus 
the dynamics of the scalar field (j) is ruled by the equation 

V^0 = 3a J^. (3) 

In vacuum (i.e., T'^'^ = J'^ = 0), or in the absence of diffusion (cr = 0), (/> is constant 
throughout spacetime and the evolution equations ([2]) reduce to the Einstein equations 
with cosmological constant. When diffusion is present, we get a "variable cosmological 
constant" model. Note however that in contrast to "the variable cosmological constant" 
theories considered in the literature], in the present situation the dynamics of the 
cosmological term is not prescribed a priori, but it is determined by the matter fields 
through the equation ([3]). 

The general relativistic theory of diffusion briefiy outlined above is developed in 
more details in sections [3]|4] and it is applied in section [5] to the cosmological model 
with the simplest geometry, namely the fiat Robertson- Walker spacetime. This model 
is simple enough to avoid serious mathematical difficulties and at the same time it 
is sufficient to derive some interesting physical implications of the new theory. In 
particular, it will be shown that the model predicts future collapse of the universe into 
a singularity in finite time even for positive (but sufficiently small) initial values (po 
of the cosmological scalar field. Since 00 can be identified with the current observed 
value of the cosmological constant A, this result is in contrast with the conclusions 
of the standard (diffusion-free) cosmological model with A > 0, which predicts an 
unlimited phase of accelerated expansion of the universe in the future. The reason for 
the formation of a singularity in the presence of diffusion is that the cosmological scalar 
field (p, although initially positive, may become negative at later times, after which the 
cosmological model behaves qualitatively like a solution of the Einstein equations with 
negative cosmological constant. 

Before going to the main topic of the paper, I shall review in the next section some 
aspects of the kinetic theory of diffusion in the non-relativistic and special relativistic 
case. 
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2 Kinetic diffusion on flat spacetimes 

Let f{t,x,p) > denote the distribution function in phase space of an ensemble of unit 
mass particles. The integral of /(t, x,p) in a region V xU of phase space is the number 
of particles which, at the time t >0, have position x G V and momentum p £ U. The 
total number of particles in the system is then 

nit) = f{t,x,p)dpdx. 

The kinetic evolution of the particles system is governed by a partial differential equa- 
tion on /; the particular form of this equation depends on the interaction among the 
particles. The most basic example is the non-relativistic free-transport equation: 

dtf+p-V,f = 0. (4) 

The free-transport equation describes the kinetic evolution of free-moving particles. If 
the particles are interacting either by internal or external forces, a new term is to be 
added in the right hand side of Q. For instance 

dtf+p-VJ = aApf, (5) 

is the kinetic Fokker-Planck (or Kramers) equation in the absence of friction |27j . 
Here a > is the diffusion constant and Ap denotes the Laplace operator in the p 
variable. The new term Apf is the most common and basic way to model diffusion in 
the momentum variable. The main physical assumption behind equation ([5]) is that 
the particles are moving in a background fluid in thermal equilibrium (thermal bath). 
Assuming that the molecules of the fluid are much lighter than the particles, and that 
the particles make up a sufficiently dilute system, the total force acting on the particles 
can be macroscopically approximated by two dominant contributions: diffusion, which 
is due to thermal fluctuations and is therefore associated — at the microscopic level — to 
random collisions with the molecules of the fluid, and friction, which takes into account 
deterministic grazing collisions among the particles. In this paper the contribution due 
to friction is neglected. 

Note that n = const., i.e., the total number of particles is conserved. This is a 
trivial consequence of the divergence theorem: 
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dtfdpdx = I I {—p ■ Vxf + o-Apf)dpdx = 0. 



Furthermore the average kinetic energy, defined as 



£km{t) = 2 I \pffdpdx, 



increases linearly in timqj, since 



3-^kin(i) = 3o-n. (6) 

dt 

*In particular, the particle system is unable to relax to an equilibrium state under the action of the 
diffusion forces alone. To this purpose, a friction term is to be added to the right hand side of ([5]), see [27] . 



The relativistic, Lorentz invariant generalization of ([5]) iqf|: 

p^5,M/ = aAW/, (7) 

where p^ — the particles 4-momentuni — satisfies the mass shell condition 

where r] is the Minkowski metric. The previous equation can be used to express the 
time component of the 4-momentum in terms of its spatial components: 



The operator in the left hand side of d?]) is the relativistic free-transport operator. 
In the right hand side, Ap denotes the Laplace-Beltrami operator associated to the 
hyperbolic metric h: 

ij — ij PiPj ) 

where p = p/p^ is the relativistic velocity. Note that h is the Riemannian metric in- 
duced by the Minkowski metric over the hyperboloid p*^ = y^l + |pp. For a justification 
of ([7]) see [1], where the more general case with friction is discussed. The explicit form 

of Ai'^V is 

aW/ = -^d,. {Vd^ih-'fd,.f) = I9,. (^^^^9,./) , 

where det /i = 1 + |pp, and (h~^y^ = 5^^ + p'^p^ is the inverse matrix of hij. It follows 
that d?]) can be written in the divergence form 

by which it is clear that n = const. Moreover, defining the average kinetic energy of 
the particles as 

Skinit) = / \/l + \p\'^fdpdx, 



it is easy to check that equation ([6]) holds in the relativistic case as well. 

The relativistic current density vector and energy-momentum tensor are given by 



J'^it,x)= [ f{t,x,p)p^%, 

7r3 p^ 

T^%t,x)= f fit,x,p)pf^p''^, 

7r3 p^ 



"■Throughout the paper, greek indexes run from to 3, latin indexes from 1 to 3 and the Einstein 
summation rule applies. Moreover, physical units are fixed such that SttG = c = 1, where G is Newton's 
gravitational constant and c is the speed of light. 



independently of the equation satisfied by the distribution function /. A straightfor- 
ward calculation using ([8]) leads to the identities 

9^pJ/^ = 0, ^^P^T^"' = ?,ar, (9) 

which imply 

a^M^^.T^^ = 0. (10) 

To conclude this section, let us comment on the fact that the particles kinetic energy 
is not conserved in the models considered above. The physical explanation for this fact 
is that, while the particles are colliding with the molecules of the surrounding fluid, the 
change of energy of the fluid molecules is neglected in the diffusion approximation, since 
the fluid is assumed to be in thermal equilibrium. If the dynamics of the coupled system 
particles-fluid were considered, then the total energy would be conserved. However, as 
shown in this section, it is possible to consider an approximation of the full dynamics 
by looking at the evolution of the particle distribution / only and still end up with 
a consistent equation for /, which is Galilean invariant in the non-relativistic case 
(eq. ([5])) and Lorentz invariant in the special relativistic case (eq. ([7])). We shall see 
that in general relativity this approximation is inconsistent. 

3 Kinetic diffusion on curved spacetimes 

Suppose now that the background Minkowski spacetime is replaced by a general Lorentzian, 
time-oriented manifold {M,g). For the moment the metric g is assumed to be given. 
Let X denote an arbitrary point of M and x^ a (local) system of coordinates on an 
open set U C M, x G U, with x^ = t being timelike. The vectors dxt^ form a basis of 
the tangent space TxM and the components of p G TxM in this basis will be denoted 
by p^. {x^,p'^) provides a system of coordinates on TU C TA4, where TM denotes the 
tangent bundle of M. The (future) mass-shell is the 7-dimensional submanifold of the 
tangent bundle defined as 

UM = {{x,p) G TM : g{x){p,p) = —1, p future directed}. 

On the subset IIL'^ = {{x,p) G UM : x G U} of the mass-shell, the condition 
g{x){p,p) = —1 is equivalent to g^uP^p'^ = —1 (where g^u = g^i^{x")), which can 
be used to express p^ in terms of p^,p'^,p^, precisely: 
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gojP^ + yiaojp'y - 500(1 + gijp'p)) 



where the choice of the positive root reflects the condition that p is future directed. 
Differentiating the mass shell relation gfj,uP^P^ = — 1 one obtains the useful relations 

^-^° = -^rV, 5..P° = -^, (11) 



which will be used below to derive some important identities. In the previous equations, 
^ fiu denote the Christoffel symbo' 
the matrix g^^ and its inverse g^" . 



^ fiu denote the Christoffel symbols of g and the indexes are lowered and raised with 



Let L denote the geodesic flow vector field on the tangent bundle. In the vector 
fields basis (S^^m , dpv ) it is given by 

The Liouville, or Vlasov, operator L is defined as the projection of L on the mass shell. 
Using (llip it is easy to derive the following local coordinates representation of L: 

The first fundamental change compared to the flat case is the definition of the "free- 
transport operator". It is now assumed, in agreement with the equivalence principle, 
that in the absence of any interaction other than gravity, i.e., in free- falling motion, the 
particles move along the geodesies of {M,g). So now the free-transport equation (j4]) is 
replaced by the Vlasov equation: 

L/ = 0, i.e., p'^5,M/-r%pV9p./ = 0, (12) 

where / : IIM — t- [0, c«) is the distribution function of particles, which is a smooth 
function on the mass-shell. See [21 \T5\ [26] for an introduction to kinetic theory and the 
Vlasov equation in general relativity. Physical systems which are supposed to be well- 
modeled by the Vlasov equation include galaxies, in which stars play the role of the 
particles, or even clusters of galaxies, where the galaxies themselves are identified with 
the Vlasov particles [7]. Kinetic theory and the Vlasov equation also have important 
applications in cosmology [8j. 

To transform (112^ into a Fokker-Planck equation we need to add a diffusion operator 
on the right hand side. Let HxM denote the fiber over x £ M of the mass-shell and 
TTx '■ HM — )■ HxM the canonical projection onto it. The action of the diffusion operator 
on / will be defined as a differential operator acting on f o tt^. The quadratic form 
g{x) induces a Riemannian metric h{x) on IIj^M. The components of h are functions 
of the coordinates {x^,p'^), which, by (llip . are given by 



Pi Pj , PiPj 



hij = Qij - -^goj - -^goi + Oooj-^- (13) 



In analogy with the diffusion operator on the mass-shell of Minkowski space defined in 
section [21 we now define the action of the diffusion operator on / by the formula 

Pp/ = AW(/ovr,), 

where Ap is the Laplace-Beltrami operator of the Riemannian metric h on HxM. The 
subscript p reminds that the operator Dp acts on the tangent space variables only. The 
expression of T>p in local coordinates is given as before by 

1 

Vdeth 



Vpf = -7==dp. (Vd^{h-^rdpj 



where now h^j = hij{x",p^) is the matrix ()13p . The generalization of ([7|) on the curved 
spacetime (M, g) is then given by 

p^dx.f - T^.p^p'^dp.f = aVpf. (14) 



The current density J^ and the energy-momentum tensor T^^^ are defined as 

J^(x) = ^\f f ^dp'^\ (15) 

T^^(x) = v^/ f?^dp^''\ (16) 



in^u -Po 

In the previous definitions, \g\ = \ det g\ and dp^'^^ denotes the 1-form dp^ A dp^ A dp^. 
Since —\/\g\/po = Vdet h, T^'^ and J*^ are, as they should be, integrals of / on the 
mass-shell with respect to the invariant measure induced thereon. Note also that, by 
the mass shell condition Qf^yP^p^ = — 1, 

TigT = g,,,Tf''' < 0. (17) 

(Equality holds for particles with zero rest mass.) Moreover it is shown in [15] that 
J^ is a timelike vector field and that T^'^ verifies the dominant and strong energy 
conditions. 

Now denote by V the Levi Civita connection of the metric g. The analogs of the 
identities ([9]) hold: 

V^J^' = 0, (18) 

^^T'"' = 3ar. (19) 

The easiest way to establish the previous identities at the arbitrary point x £ C/ is to 
use a coordinates system x^ such that g^i, = rj^^ and dx^g^v = at x. For instance, 
for (|19p we have 

V.T^^ - (V. log ^)T^^ + ^i/VV./ ^ + ^X^/ V. (-^) ..-3 

where the symbol = means that the equality holds at x G f/ in the coordinates x^. 
Integrating by parts in the last integral we conclude that (J19p holds at x in the coordi- 
nates x'^ and therefore in any other coordinates. The proof of (|18|) is similar. Finally 
by ([I8])-([I9]) the analog of equation ^ holds 



V^V^r^"" = 0. (20) 

4 Coupling with the Einstein equations 

In this section I will address the main question of this paper, which is how to couple 
the Fokker-Planck equation to the Einstein equations 

R^iv — 7;RgtJ.u = T^y. (21) 



By the Bianchi identity V^{R^i, — ^g^yR) = 0, the energy momentum tensor must sat- 
isfy V^Tfj^iy = 0, otherwise solutions of the Einstein-matter system cannot exist. Since 
for solutions of the Fokker-Planck equation (J14p the energy-momentum tensor ()16p is 
not divergence free, it is then meaningless to use the tensor ()16p in the right hand side 
of (|2ip . This incompatibility can be solved by assuming that there exist other matter 
fields in spacetimqj. The role of these additional matter fields is similar to that of the 
thermal bath in the theory outlined in section [2j Let T^i, denote the energy- momentum 
tensor of the additional matter fields and consider the Einstein equations 

Now the Bianchi identities imply, using (J19p . 

V^Jf"" = -3ar. (23) 

As a way of example, consider a system of particles undergoing diffusion in a perfect 
fluid. In this case the tensor T^j^ is given by 

'^^lu = {p + p)Ufj,Uu+pg^,u, 

where p is the rest frame energy density, p the pressure and u^ the 4-velocity of the 
fluid. Projecting (I23p in the direction of u^ and onto the plane orthogonal to u^ we 
obtain 

Vf,{pu^') + pVf.u" = SaJ^'u^,, (24) 

{p + p)u''V^Uu + UuU^'V^p + VuP = -3a{Ju + Uu{Jf,uf')). (25) 

For (T = (i.e., when the particles and the fluid are not interacting), the previous 
equations reduce, respectively, to the continuity equation and the Euler equation of a 
perfect fluid in general relativity, see f29]. 

An alternative solution to the incompatibility of equations (J14p and (j2ip is to modify 
the left hand side of the Einstein equations by adding to it a new tensor term %fj,u- 
Although this solution is formally equivalent to the previous one (by setting T^j, = 
—JCfiu), it is important to distinguish their physical meaning. The simplest choice for 
%fj_iy is the cosmological scalar field term 

^tiiy = Hfiu, (26) 

where (/> is a scalar fielqj. This choice leads to the following modiflcation of ()2ip : 

Rfii/ — -^Rdfiv + <t^gfj.v = T^y. 
The evolution equation for (j) resulting from (|19p is 

V^0 = 2,aJ^,. (27) 



^The same hypothesis is introduced in the so-called gravitational aether theories, see [23] and the refer- 
ences therein. 

^In quantum fields theory it is customary to interpret (j) as the vacuum energy. 



In the absence of diffusion, (j) is constant. In this case the model reduces to the Einstein- 
Vlasov system with cosmological constant. 

In the rest of the paper we shall focus on the cosmological scalar field model, since 
it provides a simple and yet physically interesting framework to study the effects of 
diffusion on the dynamics of spacetime and matter. 

By (I27p . diffusion in a scalar field is only possible when J^ is irrotational: 

y^.Ju - ^vJ^. = 0. (28) 

The previous equation imposes a rather severe restriction on the solution / of the 
Fokker-Planck equation (I14p . which makes the question of existence of solutions a 
highly non-trivial task. An important example where (j28p holds is the class of spatially 
homogeneous and isotropic spacetimes, which are the most popular models for the 
universe in cosmology. Cosmological models in this class are studied in the next section. 
By defining the (average) 4-velocity of the matter as 



and using ([27|) we obtain 



u^'V^cj) = -3a(-J^J^)i/2. (29) 

It follows by (j29p that the cosmological scalar field (p is decreasing along the matter 
flow. Hence the model under discussion predicts that the particles gain energy by 
diffusion and that this energy is provided by the cosmological scalar field (p. Moreover 
it follows by (pUj) that the cosmological scalar field (p satisfies the homogeneous wave 
equation: 

D,/. = 0, n = V'V^. (30) 

This implies in particular that the energy-momentum tensor of (/), 

V = V^cAV^c/. - ^9m.V",^V„(/. (31) 

is divergence free: 

V"5^, = 0. (32) 

Therefore the cosmological scalar field propagates throughout spacetime in form of 
waves without dissipation. 

5 The spatially homogeneous and flat isotropic 
cosmological model 

The purpose of this section is to study the qualitative behavior of spatially homogeneous 
and isotropic solutions of the system 

Rfiv — ■::R9fi.iy + (pdfj.!^ = Ti_i,y, (33) 

P^d,, f - r^^p^p^dp. f = aVpf, (34) 

10 



where T^^ is given by (|T6|) . In this apphcation, the particles undergoing diffusion 
represent the galaxies. 

The symmetry assumption means two things. Firstly that there exist six Killing 
vectors fields Xa, j4 = 1, . . . 6 on spacetime, whose orbits are spacelike hyper surf aces. 
Secondly that the distribution function / satisfies Xjif = 0, where Xa, A = 1, ... 6 
are the lifts over the mass shell of the Killing vectors Xa- Under the given symmetry 
assumptions, the metric can be written in the standard form 



ds^ 



-dt^ + a{tf 



dr"^ 



1 — kr^ 



+ r'^{de^ + sm'^ 



where fc = or ±1 is the curvature parameter of the hypersurfaces t = const. As 
proved in [25], the symmetry assumption on the distribution function restricts / to be 
of the form / = F{t,q), where q = {pf + p\ +p|)^'^. In particular, the 4-current J'^ 
and the energy momentum tensor T^jy take the form 

J^ = (n, 0,0,0), T^, = diag(p,?,y,y), 

where the total number of particles n is given by 

47r 



n{t) 



F{t,q)q'dq, 



aity Jo 
while the energy density p and the pressure J" are given by 

4tt 



Pit) 



a{tY 



q^F{t,q)^a{tY + qHq, 



(35) 



nt) 



A-K 



q^F{t,q) 



Observe the inequalities 



3a(t)4 Jo ^a{tf + g2 ^• 



p{t) > nit), ?{t) < ^p{t), 



(36) 



the second one being equivalent to ()17p under the assumptions of spatial homogeneity 
and isotropy. To write the Fokker-Planck equation it is more convenient to use the 
coordinates pi,P2,P3, rather than q. Setting F{t,q) = fit,pi,P2,P3), eq. ([M|) simplifies 

to 

'a{t)'^6ij+p,pj\ 



dtf = dp^ 



(37) 



where 



P 



pI+pI+ pI 



'1 + 



a{tf 



Note that in the Minkowski case (a = 1), ((371) reduces to ([8]) in the spatially homoge- 
neous case. 

The cosmological scalar field is a function of t only and the equation (|27p becomes 

6 = —San. 
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^W = -5^- (38) 



The local conservation of the number of particles, eq. (jlSp . implies 

a{ty 

where a subscript stands for evaluation at t = 0. Denoting wq = SnoOg > 0, the 
equation for cp becomes 

Note that the equation for (p contains only matter terms evaluated at time t = 0. This is 
a general consequence of the fact that satisfies the homogeneous wave equation (130p . 
From now on I shall restrict the discussion to the spatially flat case fc = 0. In this 
case it is possible to introduce "cartesian" coordinates such that the metric takes the 
form 

ds'^ = -df + a{tf{dx^ + dy^ + dz^) 

and ao = 1. The Einstein equations are 

/a\ a f d\ , ^ 

3 - -</> = /0, -2-- - +</) = ?. 
\a J a \a J 

It is convenient to rewrite them as 

H = U-\{p + ^y)-H\ (39) 

6 D 

H^:={^^' = \{p + cP), (40) 

where H = a/ a is the Hubble function. 

Initial data at t = for the system (I35])-(I10]) consist of Hq = H{Q), 4>o = (f>{0) 
and < -Fo(g) = F{0,q) = fo{pi,P2,P3) G C^(K^), such that po is bounded and the 
Hamiltonian constraint (|40p is satisfied at t = 0. It will be assumed that Hq > 0, 
since this is the most interesting case for the applications (the universe is currently 
expanding). It will also be assumed that /o does not vanish identically, since otherwise 
the only solution is vacuum. It follows from (j40p that for admissible initial data the 
value of 4>o is determined by the initial data Hq and /o through the Hamiltonian 
constraint 

</.o = mi - Po. (41) 

Finally, (pQ will be identified with the present observed value of the cosmological con- 
stant. In particular, the case (/fo ^ is the most interesting for the applications. Since 
it is an ubiquitous quantity in cosmology, I shall also use 



"» = 5% = ^-M' <^^> 



together with the following additional Hubble-normalized constants: 

Using (jiOj) . eq. (p9|) becomes H = —{p + '?)/2; in particular H is always decreasing 
and the following two possibilities may occur: 

12 



1) Either H remains positive for all times, or 

2) H vanishes at some time t* > 0. 

In the first case the universe is expanding forever in the future, in the second case it 
collapses into a singularity in finite time. Which of the two possibilities occurs depends 
on the initial value of the cosmological scalar field (p. If 00 > is such that 4> remains 
always positive, then, by (j40p . H remains positive as well. On the other hand, if 
vanishes at some time t, then by ()39p . H < — |i?i>|/3 — H^, for t > i and therefore H will 
become zero at some time t* > t. After this moment, the Hubble function continues to 
decrease until it diverges to — oo at some finite time t*^, > t^, and a{t) — )• as i — )• t^,*. 
Since (p is decreasing, future collapse in finite time occurs when 0o < 0. In particu- 
lar, even though the initial value of (po, i.e., the present observed cosmological constant, 
were vanishing, the model under study still predicts future collapse into a singularity 
in finite time. Even more remarkable is that the same conclusion holds if (po is positive 
and small. In fact it will now be shown that if 

< 00 < ^, (44) 

-Wo 

then (f> becomes zero in finite time. To see this, suppose on the contrary that (p > for 
all times. Since H is decreasing we have 

a{t) = Ha{t) < Hoa{t). 

Gronwall's inequalit}o now gives 

a{t) < e^o*. (45) 

Combining (|l5|) with ([38]) we get 



,, , , f~^°° ds , TT-ocr 

hm cp(t) = (po- ujoa / ——r < (po - -yj-, 

which is negative under assumption (I44p . Whence (p must vanish at some time i £ 
(0, +oo), a contradiction to the hypothesis that (p was always positive and the claim is 
proved. Note also that by (jH]), and using the Hubble normalized variables (ji2]) - (H3]) . 
the condition found for future collapse in finite time (including cpo < 0) can be rewritten 
as 

fio > 1 - No^o, 

which is satisfied by all admissible initial data when the diffusion constant is sufficiently 
large: 

So > Nq\ 

For So = 0, the model reduces to the Einstein- Vlasov system with cosmological con- 
stant A = (p = const and the condition for future collapse in finite time becomes fio > 1, 
which is equivalent to A < 0. The latter condition is optimal, since the Einstein- Vlasov 

*See [H]. 
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system with non-negative cosniological constant in a flat Robertson- Walker geometry 
is forever expanding in the future for all initial data. 

Next I will show the existence of initial data for the cosmological scalar field for 
which (p, and therefore H, remains positive for all times, i.e., the universe is forever 
expanding in the future. By what has just been proved, this can only happen if 
a < SHq/uq, i.e., Sq < N^ . The following argument requires a stronger bound on a, 
precisely 

a<^, i.e., So<(3iVo)-\ (46) 

no 

which will be therefore assumed from now on. Thanks to (1461) we can choose 



</'o>^, i.e., r?o < 1 - 3iVoSo. (47) 

-Wo 

By continuity, there exists < t < +oo such that (j){t) > on [0, t) and if t were finite 
we would have 0(i) = 0. Our goal is to prove that i = -|-oo. Using the bound "P < p/3 
and that H is decreasing on [0, f), we obtain 

6 b 6 6 



2 



2ij2 > _2H^. 



It follows that 



H{t) > Y^^ => a{t) > ao(l + 2Hot)l. 



for all t G [0,i). Then from (|38p we obtain 

, , . , f* ds , SnoCT 

Jo a{sy Ho 

Assuming t < 00, assumption (j47p entails (j){t) > 0, a contradiction, and therefore 
t = -|-oo must hold. Moreover, if the strict inequality holds in (j47p our model behaves 
like a perfect fluid Robertson- Walker model with negative spatial curvature {k = —1), 
since H never vanishes, not even asymptotically. 
To summarize: 

• Future collapse in finite time takes place for all admissibile initial data when 
So > A^o~^ and for l^o > 1 - A^o^o when So < N^^. 

• Unlimited future expansion takes place for initial data such that So < (3Ao)~^ 
and r^o < 1 — 3AoSo. 

More accurate estimates than those presented in this section are required to cover the 
full range of initial data. Numerical simulations may also provide a useful insight to 
the problem. 



14 



6 Summary and final remarks 

In this paper a new model to describe the dynamics of particles undergoing velocity 
diffusion in general relativity has been proposed. The main goal was to understand 
the influence of the diffusion forces on the large scale dynamics of the galaxies in the 
universe, in the approximation where each galaxy can be represented by a particle of 
the system. 

The dynamics of the particle system was described through the kinetic distribution 
function /, which was assumed to satisfy a Fokker-Planck equation without friction 
on the tangent bundle of spacetime. In the absence of diffusion, the Fokker-Planck 
equation reduces to the Vlasov equation in general relativity. 

Since the energy momentum tensor associated to a solution of the Fokker-Planck 
equation is not divergence free, the coupling with the Einstein equations is inconsistent. 
Two solutions to this problem have been proposed: Either assume the existence of 
additional matter fields in spacetime, or modify the Einstein equations. This paper 
explored in some details the case in which a cosmological scalar field is added to 
the left hand side of the Einstein equations, providing a natural generalization of the 
cosmological constant theories. It has been shown that diffusion changes drastically 
the qualitative properties of even the simplest cosmological models, namely the flat 
Robertson- Walker spacetimes. 

The paper cannot be concluded without a few comments on the physical limita- 
tions of the models presented here. First of all it should be emphasized (once again) 
that diffusion is not a fundamental interaction, but rather an approximation for the 
dynamics of two particle systems, one of which is dynamically dominant on the other. 
While in the Newtonian and special relativistic diffusion theories the dynamics of the 
"dominant component" can be broken apart from the dynamics of the other "weak 
component", the latter being therefore relegated to a background thermal bath, in 
general relativity this approximation is consistent and some interaction between the 
two components must always be taken into account. The next best "diffusion approx- 
imation" consists in assuming the simplest possible dynamics for the thermal bath, as 
the cosmological scalar field theory developed in this paper. However it seems that 
the only theoretical way to justify a specific model for the thermal bath is by deriving 
the diffusion approximation model from a (still unknown) fundamental theory for the 
interaction of the two particle systems. 

A final comment on the diffusion models introduced in this paper is that they lack 
of an action principle formulation, which is again due to their phenomenological, rather 
than fundamental, nature. 
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